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ABSTRACT 



Context, MRI turbulence is a leading mechanism for the generation of an efficient turbulent transport of angular momentum in an 
accretion disk through a turbulent viscosity effect. It is believed that the same process could also transport large-scale magnetic fields 
in disks, reshaping the magnetic structures in these objects. This process, known as turbulent resistivity, has been suggested and used 
in several accretion-ejection models and simulations to produce jets. Still, the efficiency of MRI-driven turbulence to transport large- 
scale magnetic fields is largely unknown. 

Aims. We present new analytical and numerical results aiming at quantifying the turbulent resistivity produced by MRI-driven turbu- 
lence in accretion disks. 

Methods. We investigate this problem both analytically and numerically. We introduce a linear calculation of the MRI in the presence 
of a spatially inhomogeneous mean magnetic field. We show that, in this configuration, MRI modes lead to an efficient magnetic field 
transport, on the order of the angular momentum transport. We next use fully non linear simulations of MRI turbulence to compute 
the turbulent resistivity in several magnetic configurations. 

Results. We find that the turbulent resistivity is on the order of the turbulent viscosity in all our simulations, although somewhat lower. 
The variations in the turbulent resistivity are correlated with the variation in the turbulent viscosity as a function of the imposed mean 
field. Finally, the turbulent resistivity tensor is found to be highly anisotropic with a diffusion coefficient 3 times greater in the radial 
direction than in the vertical direction. 

Conclusions. These results support the possibility of driving jets from turbulent disks; the resulting jets may not be steady. 
Key words, accretion disk - turbulence - mhd 



1. Introduction 

After several decades of debate, and in spite of significant 
progress, the origin and efficiency of angular momentum trans- 
port remains a central problem in accretion disks physics. 
Turbulence is one of the major physical mechanism through 
which the observationally well-constrained anomalous transport 
can b e achieved; indeed, the first a model (Shakura & Sunvaevl 
1 1973b already assumed a strong level of turbulence, leading to 
an effective viscosity orders of magnitude higher than molecu- 
lar viscosity. However, the physical origin of this turbulence in 
disks is still largely discussed. 

Various sources of turbulence, both hydrodynamical and 
magneto- hydrodynamical, have b een proposed. Subcritical tur- 
bulence ( Richard & Zahnl 1 19991 and references the r ein), if 
prese nt, is too inefficient dLesur & LongarettH 120051 ; iJi et al.l 
2006). Convection is also too inef ficient and transports angular 
mom entum in the wrong direction dCabotHl996HStone & Balbusl 
119961) . Two-dimensional turbulence driving by smal l-scale, in- 
coherent gravitatio nal instabilities jGammid [T 996) or baro- 
clinic instabilities dKlahr & Bodenheimert 120031) is a possi- 
ble option, although the last one is still high ly controversial 
dJohnson & Gammidl2006l;lPetersen et alj|2007l) . 

In a seminal paper, iBalbus & Hawlevl d 199 lb indentify an 
MHD instability, the magnetorotational instability (MRI) that 
drives turbulence in the nonlinear regime. This instability pro- 
vides the most extensively studied transport mechanism, mainly 



with t he help of local un stratified dHawlev et al.ll 19951) a nd strat- 
ified dStone et al.l fl996l) 3D simulations, and global dHawlevI 
2000) disk simulations. These simulations have shown that MRI 
turbulence is an efficient way to transport angular momen- 
tum, although the rol e of microphysical processes has largely 
been underestimated dLesur & Lo ngaretti 2 0071: iFromang et al.l 
l2007h . 

MRI turbulence may also produce resistive transport (trans- 
port of magnetic fields through a "turbulent resistivity" process) 
in disks. This tra nsport is a key ingredient of accretion- ejection 
models (see, e.g.. lFerreiralll997l;ICasse & Ferreirall2000l and ref- 
erences therein). The related turbulent resistivity r\j is parame- 
terized with the Shakura-Sunyaev ansatz as r\j - u^VaH (Va is 
the Alfven speed based on the mean field); stationary accretion 
ejection models require an anisotropy of the turbulent resistivity 
transport of few; furthermore, a very efficient turbulent trans- 
port with a n < 1 is required to produce stationary structures. 
Less efficient transport would not prevent the launching of jets 
from disks, but these structures would then be nonstationary. The 
question of whether turbulent resistivity is an adequate source of 
field line slippage through the gas to launch a jet is still contro- 
versial in the literature, and is further discussed in the concluding 
section of this paper. 



1 In this paper, a transport process refers to any physical mechanism 
by which a magnetic field line is moved from one place to another, 
without assuming anything about flux-freezing conditions. 
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In any case, turbulent resistivity is an issue in its own right, 
and in this paper, we present new numerical results aimed at 
quantifying more precisely the resistive transport due to MRI 
turbulence. Our approach differs from the method recently intro- 
duced by Brandenburg, Radler, Schrinner and coworkers to com- 
pute the turbulent resist ivity tensor ("see ISchrinner et ail 12005, 
Brandenburg et al. 2008 and references therein). Their method 
relies on an expansion of the mean electromotive force with re- 
spect to the mean field and its derivatives, and makes use of ap- 
propriately chosen test field, whose evolution is computed along 
with the tur bulent flow. Our approac h is also different from the 
one used bv lGuan & Gam mie (2009), in which a magnetic struc- 
ture is imposed as an initial condition, and the resistive proper- 
ties of the background turbulent flow are deduced from the decay 
time of this structure. The method we adopt here is described in 
section 2, and the physical content of the two methods is dis- 
cussed there. More generally, section 2 describes the physics and 
the numerical methods we have used to study turbulent resistiv- 
ity in disks. Section 3 presents a linear analysis of the model, 
which sheds some light on the numerical findings described in 
section 4. Section 5 discusses the astrophysical implications of 
these results along with some future line of work. 

2. Shearing-box equations and numerical method 

MRI-related turbulence has been extensively studied in the lit- 
erature, but little attention has been devoted to the question of 
turbulent-driven resistivity yet. One of the major advances of 
the recent years has been the realization that an accurate de- 
termination of turbulent transport properties requires an accu- 
rate representation of all scales down to the dis sipation scales 
(iLesur & Longarettill2007t iFromang et"ai1l2007l) . although this 
accuracy has not yet been achieved for astrophysically relevant 
magnetic Prandtl and Reynolds numbers regimes. We do not 
imply that actual astrophysical Reynolds numbers need to be 
resolved in simulations to obtain reliable results, an obviously 
hopeless task anyway; but for large enough Reynolds and mag- 
netic Reynolds numbers, one expects the dissipation scale to be 
decoupled from the transport scales, in which case one could 
use a closure model without resolving the dissipation scales. 
However, today simulations don't seem to have reached this 
regime yet. 

Shearing box simulations offer a particularly convenient set- 
ting to quantify turbulent resistivity; however, the boundary con- 
ditions prevent the existence of the large scale gradients of mean 
magnetic field required for characterizing the transport of this 
mean field produced by turbulence. We bypass this difficulty by 
a prescription whose physical motivation and formulation will 
be described later on in this section. For the time being, we 
briefly recall the basic equations for the shearing-box model, 
which has been larg ely studied an d used i n the lit e rature . The 
reader may consul t lHawlev et al.l d 1995b . iBalbusI d2003l) and 
Regev & Umurhaij d2008l) for an extensive discussion of the 
properties and limitations of this model. 

Since MHD turbulence in disks is essentially subsonic, we 
will work in the incompressible approximation, which allows 
us to eliminate sound waves and density waves from the prob- 
lem. Alt hough density waves are excited in shearing box tur- 
bulence (Hein emann & Papaloizou 2008b) they should not have 
a big impact on the turbulence spectrum itself or on turbulent 
transport as their am plitude decays exponentially fast when one 
goes to small scales dHeinemann & Papaloizoul2008al) . This has 
also been confirmed by direct compari son between incomp ress- 
ible and compressible simulations (Froma ng et al.l 12007). We 



also neglect vertical stratification, consistent ly with the local 
shearing-box model dRegev & Umurhan 2008). Explicit molec- 
ular viscosity and resistivity are included in our description. 

The shearing box equations follow from a local approxima- 
tion. We chose a Cartesian box centered at r — Rq, rotating with 
the disk at angular velocity ft = Cl(Ro) and having dimensions 
(L x , L y , L z ) with L, <k Rq. Assuming R^ — > x and r - Rq — > —y, 



one eventually obtains the shearing box equations: 

d f U + V • (17 ® U) = - Vn + V ■ (B <g) B) 

-2ft X U + 2Q.Sye y + vAU, (1) 

d,B = V X(U X B) + i]AB, (2) 

V • U = 0, (3) 

V • B = 0. (4) 



The boundary conditions associated with this system are peri- 
odic i n the x and z directi on and shearing-periodic in the y direc- 
tion dHawlev et alj fl995). In these equations, we have defined 
the mean shear S = -rd r Q., which is set to S — (3/2)Q, assum- 
ing a Keplerian rotation profile. The generalized pressure term 
II includes both the gas pressure term P/po and the magnetic 
one B 2 /2po. This generalized pressure II is actually a Lagrange 
multiplier enforcing equation (01, and is therefore computed by 
solving a Poisson equation. Note also that the magnetic field is 
expressed in Alfven-speed units, for simplicity. 

The steady-state solution to these equations is the local 
Keplerian profile U = Sye x . In this paper, we will consider only 
the turbulent deviations from this Keplerian profile. These may 
be written as V = U -Sye x , leading to the following equations 
forV: 



d,V + V • (V <g> V) = - Vn + V • (B ® B) - Syd x V 

+(2ft - S)V y e x - 2QV x e y + vAV,(5) 



d,B = -Syd x B + SB y e x 

+V X (V X B) + t]AB, (6) 

V • V = 0, (7) 

V • B = 0. (8) 



Following lHawlev et al.l d 1995b . one can integrate the induc- 
tion equation (O over the volume of the box, leading to: 

^=S(B y )e x , (9) 
ot 

where () denotes a volume average. Therefore, the mean mag- 
netic field is conserved, provided that no mean radial field is 
present. In this work, the mean field will be either vertical or 
azimuthal. 

To numerically solve the shearing-box equations, we use 
a spectral Galerkin representation of equatio ns ([I])-® in the 
sheared frame (see ILesur & Longarettill2005l) . This frame al- 
lows us to use a Fourier decomposition since the shearing- 
sheet boundary conditions are transformed into perfectly pe- 
riodic boundary conditions. Moreover, this decomposition al- 
lows us to conserve magnetic flux to machine precision with- 
out any modification, which is an advantage compared to finite- 
difference or finite-volume methods (the total magnetic flux cre- 
ated during one simulation is typically 10~ n ). Equations © and 
([8j a re enforced to machine precision using a spectral projec- 
tion (Pevrel l2002l) . The nonlinear terms are computed with a 
pseudospectral method, and aliasing is prevented using the 3/2 
rule. To always compute the physically relevant scales in the 
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sheared fr ame, we use a remap met hod similar to the one de- 
scribed bv lUmurhan & Regevl d2004l) . This routine redefines the 
sheared frame every T Kma p = L X /(L V S) and we have checked 
that none of the behaviour we describe in this paper was related 
to this numerical timescale. Since spectral methods are very lit- 
tle dissipative by nature, we check that numerical dissipation is 
kept to very small v alues, computing the total energy budget at 
each time step (see iLesur & Longaretti 2005, for a discussion 
of this procedure). We therefore ensure that numerical dissipa- 
tion is responsible for less than 3% of the total dissipative losses 
occurring in these simulations. 

To quantify the dissipation processes in the simulations, we 
use dimensionless numbers defined as: 

- The Reynolds number, Re — S L%/v, comparing the nonlinear 
advection term to the viscous dissipation. 

- The magnetic Reynolds number, Rm = SL^/rj, comparing 
magnetic field advection to the Ohmic resistivity. 

- The magnetic Prandtl number, defined as the ratio of the two 
previous quantities Pm = Rm/Re = v/77, which measures 
the relative importance of the dissipation processes, and, cor- 
relatively, is related to the ratio of the viscous and resistive 
dissipation scales. 

In the following we use S~ l as the unit of time and SL Z as the 
unit of velocity. One orbit corresponds to T m \, = 3nS~ l . For 
simplicity, we keep the same notation for dimensionless and di- 
mensional quantities. 

2.1. Turbulent resistivity definition 

Our aim is to test to which extent the effect of MRI turbulence 
on the mean field can be modelled as a turbulent resistivity on 
large scales, and characterize inasmuch as possible the result- 
ing turbulent resistivity tensor. We therefore distinguish the large 
scale mean field B and velocity V and the fluctuating (turbu- 
lent) fields b and v. We assume (b) = and (v) = where 
() denotes an ensemble (or time, under the ergodic hypothesis) 
average. The induction equation reads: 

d,B~ = V X (Vxfi)+Vx£ (10) 

where we have defined the mean electromotive force (EMF) £ = 
(v x b). The turbulent resistivity hypothesis assumes: 

8i=/3 m d Xk Bj. (11) 

Our main objective is to test and quantify this assumption. We 
keep no term proportional to B (a-like effect), as we found out 
that they are not produced in our simulations. 

At this point, two different routes are open to study the tur- 
bulent diffusion of the magnetic field. 

- In the first one, Eq. ( fTTT i is assumed. One can split Eq. (O into 
an equation for the mean field B and one for the deviation 
from the mean b. As the induction equation is linear in the 
field, b formally depends on the mean and turbulent veloc- 
ity fields and mean magnetic field only. One then introduces 
extraneous magnetic fields V q , whose equation of evolution 
is the same as for 6, except for the mean field B which is re- 
placed by a adhoc test field B 11 . It is then possible to derive 
the EMF associated to this test field G pq = (v x b pq ). 

One can deduce some of the components of B,^ using several 
test fields and computing the correlations between B pq and 



S pq . By construction, this method probes the velocity field v 
with a test field B pq which is different from the real field B 
entering in the Lorentz force. Therefore, one assumes implic- 
itly that the properties of /3y* don't depend on the topological 
properties of B , which might not be tr ue for subcritical dy- 
namos (see e.g. lLesur & Ogiiviell2008l) or turbulence driven 
through the L orentz force like MRI turbu lence. This method 
is deta iled in ISchrinner et alj d2005l) and iBrandenburg et al] 
(2008) and references therein. 
- Alternatively, one can impose the constancy of some large 
scale component of the magnetic field gradient throughout 
the evolution and test if Eq. ( fTTt accurately represents the 
effect of turbulence on field diffusion. As in the previous ap- 
proach, adequate choices of the imposed field configuration 
allow us to characterize various elements of the diffusivity 
tensor. This approach is adopted here and specified in more 
detail in the following subsection. 

These two approaches do not quite probe the same physical 
aspects of the problem. This first one is justified inasmuch as 
the mean field hypothesis Eq. (TTTT) holds, provided that the tur- 
bulent resistivity tensor does no depend on the field topology. 
The second one is justified without assumption if the turbulent 
diffusion of the field is much slower than the other time-scales 
associated to the turbulence dynamics. However, even if this re- 
quirement is not satisfied, the physical content of this second 
approach can be grasped by comparing the questions of mag- 
netic field and momentum transfer. Consider for example a sim- 
ple unmagnetized plane Couette flow. With shearing box bound- 
ary conditions, the turbulent flow is unable to react of the shear 
(velocity gradient) impos ed on the flow through the boundary 
conditions. As shown by ILesur & Longarettil j2005), the flow 
nevertheless becomes turbulent above some transition Reynolds 
number R c and produces a shearwise transport characterized by 
R c . Alternatively, one can let the flow react on the imposed ve- 
locity gradient by adopting e.g. wall boundary conditions in 
the shearwise direction. This changes has little effect on the 
transition Reynolds num bers (at least for cyclonic rotation: see 
ILesur & Longarettill2005l) . and does not appear to affect much 
the dim ensionless transport , as su ggested by the numerical re- 
sults of ILesur & Longarettil d2005) and the experimental results 
compiled in lDubrulle et a l. (2005). These considerations make it 
plausible that the turbulent diffusion of the magnetic field with 
an imposed field gradient or with a self-consistently evolving 
one should be similar, although a more systematic study is re- 
quired to ascertain this statement. 

In practice, we chose field configurations such that only one 
component of the field gradient is non vanishing at any given 
time, so that the current is directly proportional to this non- 
vanishing gradient. In this situation, one can dispense with the 
third order tensor and use 

Si = -rjijJj (12) 

instead, where J — V X B is the mean electric current and 77 is 
a constant tensor, but the reader should keep in mind that there 
is no reason that 77,^ should be antisymmetric in the relevant 
indices. Similarly, 77 has no reason to be diagonal for MRI tur- 
bulence. In this work, we explore this model looking for linear 
correlations between the components Gj and and computing 
the related resistivity tensor coefficients. 

To keep in line with the logic of the a parameterization 
as adapted to an incompressible fluid (see ILesur & Longarettil 
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120071) . one may introduce an a tensor relating to the turbulent rj 
tensor: 



With our choice of units (5=1 and L, = 1), one has a ViJ = rjij, 
and the distinction between the two quantities is usually dropped 
in the remainder of the paper. 

2.2. Numerical protocol 

To compute the turbulent resistivity, one needs a mean current, 
which is not present in local (shearing box) simulations. We ob- 
tain this current by imposing a large scale and non homogeneous 
field in the box. This is done with the help of our Galerkin rep- 
resentation by forcing the largest Fourier mode in one direction 
to a constant value ABq. To trigger the MRI, we impose a mean 
field Bq which can be either azimuthal or vertical. The aspect ra- 
tio is L x xL y xL z =4x2x 1 where x — <p, y = -r and z-z. The fac- 
tor 2 in L y allow us to trigger more easily secondary instabilities 
in t he strong mean ve rtical field cases (see lGoodman & Xul 19941 
and lBodo et al.l2008l) . The resol ution used is 128x 128x64 , sim- 
ilar in cell size to the one used in lLesur & Longarettil d2007). The 
Reynolds number is kept constant at Re = 1600 as well as the 
magnetic Prandtl number, which is fixed at unity: Pm = 1. This 
corresponds to a Elsasser number A = V^/r/Q. = 24 for Bo = 0.1 
and A = 1.5 for Bq = 0.025. Each simulation is integrated over 
500 shear times (to obtain meaningful time averages), and the 
average are computed from 400 last shear times (to avoid initial 
conditions artifacts). 

To postprocess the results, we first compute the time average 
of B and &j. We then use a script which extracts the mean cur- 
rent and compute the correlation with the emfs, giving in the end 
one component of the rj tensor for one run. Note that using this 
procedure, we can compute resistivity associated with B z (y) (ra- 
dially varying vertical field), B x (y) (radially varying azimuthal 
field) and B x (z) (vertically varying azimuthal field) configura- 
tions. Imposing a constant radial field B v would lead to a linear 
growth of the azimuthal field B v , which may not have any physi- 
cal meaning. We cannot compute turbulent resistivity associated 
with non axisymmetric structures because of the shearing-sheet 
boundary conditions. 

The runs we have performed do not test all the dependencies 
of the dimensionless turbulent resistive transport with respect to 
all the dimensionless parameters of the problem, but only a sub- 
set of them, namely: 

- The role of the mean field orientation Bo, and of the orienta- 
tion the imposed field varying component AZ?o- 

- The role of the relative amplitude of the mean and im- 
posed varying field components e = ABq/Bq. In the process, 
only a few components of the turbulent resistivity tensor are 
probed. 

- The role of the relative mean field amplitude, characterized 
b^B = (4npS 2 H 2 )/Bl = Bf (the last equality follows from 
our choice of units). 

No characterization with respect to the dissipation parame- 
ters has been attempted (most notably the magnetic Prandtl num- 
ber), a problem requiring substantial efforts and that will be ad- 
dressed elsewhere. 

2 Our p parameter is akin to the usual plasma beta parameter, as both 
differ only by a factor of order unity in a stratified disk in vertical hy- 
drostatic equilibrium. 
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Fig. 1. Mean emfs and best fit with the current for B = 100 and 
e = 0.3 for the fastest growing channel mode with an imposed 
vertical mean field and radially sinusoidally varying imposed az- 
imuthal current. The relevant emf component is seen to be- 
have as a resistive term. This component measures the diffu- 
sion of the vertical field in the radial direction. The amplitudes 
are arbitrary. 



Our local simulation box physical size is on the order of the 
scale height at most. As such, one would expect to mimic the ef- 
fect of large scale gradients of the magnetic field by small values 
of e. However sizeable (but of limited extent) local field gradi- 
ents may also be present in disks as a result of the dynamics, and 
it is of some interest to explore values of e of order unity as well. 

3. Analytic results 

It is of some interest to look into the linear behavior of the in- 
stability first, as it does provide some insight into the turbulent 
problem that we discuss in the next section. Indeed, there is evi- 
dence that the channel mode (both a linear and nonlinear solution 
to the incompressible equations) plays some role in the overall 
transport properties of the MRI-driven turbulent states. 

For this reason, we want to test the turbulent resistivity hy- 
pothesis (equations [10] and [TTJ in the linear regime. This is done 
assuming the fluctuations v and b are infinitely small compared 
to the averaged fields. We then calculate v and b by a classi- 
cal linear analysis and deduce the EMF £ which is a quadratic 
quantity in this limi£| . The same procedure can be used for 
Maxwell and Reynolds stresses, allowing one to compute "ef- 
fective" transport coefficients in the linear regime. 

To make comparisons easier, we investigate one of the con- 
figurations which is simulated in the next section of this paper, 
namely we impose a vertical mean field of the form 

B r (y)=B [l+ecos(A:oy)], (14) 

where Icq - 2n/L (= 2n in our choice of units) and e measures 
the relative amplitude of the sinusoidal component of the mean 
field, and is to be taken as a small expansion parameter later on 
in some of the analytic expressions. Other configurations would 
lead to a similar qualitative behavior, which is the point of inter- 
est here. 

Assuming that the development of the instability is fast com- 
pared to the viscous and resistive time-scales, we neglect the 
corresponding terms in Eqs. © and © in this linear analysis. 
Under this approximation, the system Eq. ((5]) to (0 admits a sim- 
ple equilibrium with constant pressure and varying equilibrium 

3 For this reason, this approach is also called "quasi-linear". 
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poloidal magnetic field - first order 



poloidal magnetic field - zeroth order 



Fig. 3. Zeroth and first order components of the poloidal magnetic and velocity fields for the fastest growing channel mode (B = 100, 
e = 0.3). The coupling of zeroth (resp. first) order component of the velocity field to the first (resp. zeroth) order component of the 
magnetic field gives rise to the transport of magnetic field by the channel mode. 




First order linear deviations from the mean stationary solution 
described above are referred to as XJ\, B\, P\, and one considers 
only axisymmetric perturbations of the form 



X\(y, z, t) - x(y) expO'wf - ikz), 



(17) 



where X, x refers to the velocity, magnetic, and gaz pressure 
fields. 

It is useful to distinguish the poloidal and toroidal compo- 
nents of the velocity and magnetic fields: 



Fig. 2. Ratio of the mean emf Eq and mean Reynolds stress ten- 
sor (r, (p) component due to the most unstable channel mode. 

azimuthal velocitjQ profile to balance the magnetic pressure gra- 
dient: 



V(y) = Sy 



1 dBi 



4O dy 

This leads us to introduce the total shear S t- 
1 d 2 B 2 T 



S T =S - 



4Q dy 2 



(15) 



(16) 



It turns out that only channel-like mode present correlations 
of the EMF that behave as expected under the resistivity hy- 
pothesis, so we specialize from the onset to that type of modes. 



4 An actual disk is indeed expected to react to the presence of a super- 
imposed magnetic pressure gradient at some initial time by first adjust- 
ing its velocity profile instead of its gas pressure profile. Nevertheless, 
we have also explored the converse situation, where the magnetic pres- 
sure variation is balanced by a variation of the gas pressure, without 
modification of the Keplerian profile (a behavior expected in an incom- 
pressible fluid). This change of equilibrium has little effect on the results 
reported in this section. 



v = v ± + v x e x , 
b = b ± + b x e x . 



(18) 
(19) 



Because, of the assumed axisymmetry, on can introduce po- 
tential fields \jj and <p for these poloidal components such that 



dili 

v ± = ik<j/e x H e z , 

dy 

dip 

b ± — ik<pe y H e 7 , 

dy 



(20) 
(21) 



All variables can be expressed in terms of ip, which leads to 
a second order ODE for iff. One obtains: 



(w 2 - VlKW - 2V 2 b k 3 z iff' 



(0 f _ vlfy _ 



k z t cj z + 2Q. S T V 2 k 2 



co 2 - V\k\ 



k 2 Jj = 0, 



(22) 



where VaCj) and Vt(y) are local Alfven speed type quantities, 
and Kj{y) is a local epicyclic frequency, defined by 



6 



G. Lesur and P-Y. Longaretti: Turbulent resistivity evaluation in MRI generated turbulence. 



Vl = B 2 T , 



BtB', 



Kj - 2<W()(20o - St)- 



(23) 
(24) 
(25) 



The magnetic field b and the azimuthal velocity v x can be 
expressed in terms of if/ as follows: 



b x = - 



k-B 



by = 



k 2 B T 
—i— — i/i, 



; k.SrBr 

Vx ~ tf>, 



b 7 = - 



CO 

k-Bj 



k-B' 



,,2/c in.-i _ c„t/2;,2 



u> z (S T -2n Q )-s T v z A ki 

co(co 2 - V 2 A k 2 z ) 



(26) 
(27) 
(28) 
(29) 



We have solved Eq. ( l22l numerically with periodic radial 
boundary condition^ Analytic solutions can be found in terms 
of series expansion in e in the form ifr = i//q + + ifrz • ■ ■ and 
oj 2 — u>q + co 2 + co 2 , ■ ■ .; the first order solution is rather straight- 
forward to derive, but the explicit form of this solution is not par- 
ticularly illuminating and will not be given in detail, except for 
some relevant features that will be mentioned whenever needed. 

Solutions in the e = limit are the usual MRI modes, i.e. 
simple sinusoidal modes if/o oc exp(+ik y y). Solutions for arbi- 
trary e generalize these modes. It turns out that only the node- 
less mode (k y = 0) produces a mean emf that correlates with 
the imposed current, and in the following, only these modes are 
studied. 

For definiteness, we show results pertaining to the fastest 
growing mode, for B = 100 and e = 0.3. Similar results are ob- 
tained for other sets of parameters, and for other purely vertical 
modes. 

Following the quasi-linear analysis procedure, we define the 
quadratic emf as: 



<£(y)> = i*. 



v X 



(30) 



where () denotes in this case a vertical average of the fluctu- 
ations, and the dagger stands for complex conjugation. As the 
channel mode is exponentially growing in time, the exponential 
prefactor has been left out (or equivalently, a time-dependent £ 
would be required in the above equation). The various compo- 
nents of the emf due to the unstable mode, as well as the best fit 
of Eq. ( fTTT i with the imposed current are shown on Fig.Q] 

It is of some interest to compare the efficiency of the "re- 
sistive" and "viscous" transport due to the instability, i.e., the 
value of /700 = emf I current deduced here to the value of 
v r $ — y — (Reynolds + Maxwell stress)/ shear due to the same 
mode. This is shown on Fig.|2]for the same mode and for two dif- 
ferent values of B and a range of values of e. Note that this ratio 



3 No boundary condition is implied vertically so that the modes are 
not discretized in this direction. In this section, the unit of length is set 
to Ly, but this has no incidence on the results or on the rest of the paper. 

6 The reason for this can be understood by looking at the equation 
for the first order correction ifi\, which shows that the correct radial 
behavior for coupling occurs only in this case. 



is constant throughout the growth of the mode. The analytic re- 
sult for this ratio in the limit of small e is rj^ /v r d = 27/ 1 6 - 1 .7 . 
Apparently, the relative efficiency of resistive and viscous trans- 
ports due to the instability are of similar magnitude, more or less 
independently of the amplitude of the imposed field variation 
and mean field strength. 

The comparable efficiency of the two transport effects due to 
the channel mode comes somewhat as a surprise. From a ther- 
modynamical point of view, the instability sets in because of the 
imposed velocity shear and tries to restore a thermodynamical 
equilibrium by transporting momentum to reduce this shear. At 
first sight, there is no reason that it should also try to suppress a 
gradient in magnetic field with about the same efficiency. 

However, if one considers the second-order correlation ap- 
proximation used in mean field electrodynamic, one finds (e.g. 
Radler & Rheinhardt 200^): 



Pij ~ \(v 2 )r, 



(31) 



where t is a typical correlation timescale of the small scale tur- 
bulence. This kind of result can be understood as a mixing length 
theory for magnetic fields involving velocities of the order of 
V and a length on the order of TV. Although the second-order 
correlation approximation is not entirely justified in the quasi- 
linear theory used here, at least some of the terms of the in- 
duction equation provide a mixing length behavior. This sug- 
gests that with a typical timescale on the order of the shear time, 
the dimensionless turbulent resistivity should be of the order of 
the turbulent transport (the turbulent transport being of the or- 
der of the turbulent kinetic energy). Pushing this argument one 
step further, we note that in general in numerical simulations, 
(b x by) ~ 3(v t v v ), implying that one should observe B,j ~ a/4. 
This is not consistent with the analytic result, but it agrees very 
well with the numerical results of the next section for a vertical 
field and field variation; however it does not capture the origin 
of the significant anisotropy of the resistivity tensor observed in 
these simulations. 

The analytic results as well as the preceding argument imply 
that the motion and magnetic field deviations due to the varying 
field should be coupled. An illustration of the efficiency of these 
couplings leading to this high resistive transport can be obtained 
from the analytic solution in the small e limit. In this limit, the 
correlation arises from the correlation of zeroth and first order 
quantities in e, (v x b) — (v x b\ + V\ x 6q). Fig[3]shows the 
four involved poloidal fields for the most unstable channel mode 
(B = 100, e = 0.3). Direct inspection indeed shows that the 
relevant fields do exhibit a significant level of correlation. 

Having gained this understanding of the transport due to the 
linear modes of the instability, we now turn to the full nonlinear 
problem. 



4. Numerical results 

We have simulated three different configurations, which com- 
bine a mean vertical and/or azimuthal field with a varying field 
component in the vertical or azimuthal direction, the direction of 
variation being either vertical or radial. This allows us to charac- 
terize a number of components of the turbulent resistivity tensor. 

To classify our numerical results, we use the following nam- 
ing convention: the first letter indicates the direction of the mean 
field, the second letter indicates the direction of the varying field 
and the third letter corresponds to the direction of the spatial de- 
pendency. We finally add a number to identify the simulations 



G. Lesur and P-Y. Longaretti: Turbulent resistivity evaluation in MRI generated turbulence. 



model 










a- 


ZZY1 


n 1 


0.01 


3.6 x 10" 


2 


2.9 x 10" 


2 


1.5 x 10 ' 


77Y? 


n 1 

U. 1 




3.7 x 10" 


-2 


8.5 x 10~ 


2 


i ^ v in -1 

l.J A 1U 


ZZY3 


0.1 




3.6 x 10" 


-2 


6.8 x 10" 


2 


1.5 x 10" 1 


ZZY4 


0.1 


0.04 


2.4 x 10" 


-2 


5.6 x 10- 


2 


1.3 x 10" 1 


ZZY5 


0.1 


0.08 


1 1 v ID" 

1.1 A 1U 


-2 


1 1 v 1 

1.1 A 1 U 


2 


7.8 x 10~ 2 


ZZY6 


0.025 


0.0025 


4.1 x 10 


3 


2.6 x 10- 


2 


3.7 x 10~ 2 


ZZY7 


0.025 


0.005 


1.2 x 10 


-2 


-1.2 x 10 


-2 


3.3 x lO- 2 


ZZY8 


0.025 


0.0075 


7.0 x 10 


-3 


2.8 x 10- 


2 


3.7 x 10~ 2 


ZZY9 


0.025 


0.01 


5.1 x 10" 


-3 


2.4 x 10- 


2 


3.5 x 10~ 2 


ZZY10 


0.025 


0.015 


8.4 x 10" 


-3 


1.5 x 10- 


2 


3.9 x 10" 2 


ZZY11 


0.025 


0.02 


8.5 x 10 


-3 


1.5 x 10" 


2 


3.3 x 10~ 2 



Table 1. Main results from the B~(y) (vertical field with a radial 
dependency) case in the presence of a mean vertical field 



done in the same configuration but with different field intensity. 
For instance, run XZY7 have a mean field in the x direction, a 
varying field in the z direction with a spatial dependency in the 
y direction. 

The various configurations and results are presented in the 
corresponding subsections. They are then compared in the con- 
cluding subsection. 



model 
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A 




T] zz 


a 


ZXY1 


0.1 


0.0 


0.01 


9.0 x 10~- 


1.4 x 10 J 
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U. 1 


n n 

u.u 
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1 ft v 1 1V 1 

1 .U A 1U 


ZXY3 


U. 1 
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U.UJ 


e c v 1 c\- 2 

O.O A 1U 


1.5 x 10" 1 


ZXY4 


0.1 


u.u 


0.04 


c s v 1 cr 1 

O.O A 1U 


1.4 x 10" 1 


ZXY5 


U. 1 


u.u 


0.06 


8.2 x 10~ 2 


1.6 x 10 _1 


ZXY6 


U. 1 




u.uo 


8 f) v 1 fr" 2 

O.U A 1U 


1.5 x 10 _1 


YYYI 
/\ /\ l l 


u.u 


U. 1 


U.U1 


7 A v 1 f)~ 3 

/.HA 1U 


l .U A 1U 


XXY2 


0.0 


0.1 


0.02 


1.6 x 10~ 2 


8.1 x 10 -3 


XXY3 


0.0 


0.1 


0.03 


1.4 x 10~ 2 


9.1 x 10~ 3 


XXY4 


0.0 


0.1 


0.04 


1.4 x 10~ 2 


9.2 x 10~ 3 


XXY5 


0.0 


0.1 


0.06 


1.3 x 10~ 2 


8.9 x 10^ 3 


XXY6 


0.0 


0.1 


0.08 


2.0 x 10~ 2 


1.3 x 10~ 2 


XXY7 


0.0 


0.1 


0.10 


2.1 x 10~ 2 


1.6 x 10 -3 



Table 2. Main results from the B x {y) (toroidal field with a radial 
dependency) case for several toroidal and vertical mean fields 
B°. 



The resulting ratio 77/ v (v = a with our choice of units) is 
substantially smaller than its linear counterpart, and consistent 
with the argument exposed in the previous section. 



4.1. Vertical field with radial dependency B z (y) 

In this case, we force the following structure for the mean mag- 
netic field: 



(B z ) = B° Z +6B° cos 



The mean current is then written: 



(J) = {J x )e a 



-<5fi°^sin(^)e a 



Ly 



(33) 



We are therefore looking for a correlation between J x and 
£. We show on Fig. [4] an example of a simulation result with 
B° z = 0.1 and 5B l> = 0.04 (e = 0.4). The profiles are computed 
from an average in time and in the (x,y) plane. From this figure, 
one get a classical diagonal resitivity term of j] xx ~ 5.6 x 10 -2 . 
We also find an non diagonal term r\ yx — 8 x 10~ 2 . We have 
repeated this kind of experience for various sets of parameters 
as summarized in Tab. Q] 

The extraneous EMF component was already present at the 
linear level, as noted in the previous section, but with the oppo- 
site sign. We have not been able to understand the significance of 
this sign difference. In our simulations, this component plays no 
physical role (its rotational vanishes). This may not be a generic 
feature, in particular in more complex stratified settings. This 
component shows anyway that the turbulent state is anisotropic. 

One may note a decrease of the turbulence efficiency in mod- 
els ZZY4-ZZY5, which may be due to the fact that increasing 
6Bc to high values leads to strong modification of the back- 
ground field. Since B ( ' = 0.1 corresponds to the maximum of 
the growth rate for the k z = 2njL z mode, increasing 6B° always 
weakens the instability, which might explain these results. As 
expected, this effect is not observed in B z = 0.025 models, for 
which this explanation doesn't hold. 

According to these results, and anticipating on the result of 
the last subsection where the correlations are compared in detail, 
we point out that, on average 



4.2. Toroidal field with a radial dependency B x (y) 
In this case, we impose the following field: 



(32) (B z ) = B° z 



(B x ) = B° + ( $B°co S (|^) 



The mean current is then defined by: 
(J) = W e,=^|sin(^) 



(35) 
(36) 



(37) 



As in the previous case, we look for a correlation between J z 
and £. We show on Fig. [5] an example of such a correlation for 
B° = 0.1, B° x = and 6B° = 0.02, from which we get t] zz ~ 0.08. 
Note that the other components of the EMF are poorly corre- 
lated with the imposed large scale field, indicating that either no 
correlation exists or that the turbulent viscosity model does not 
capture its physics; as a consequence, we have not quantified off- 
diagonal components of the resistivity tensor. As previously, we 
reproduce this kind of run for a range of parameters (see Tab.|2]i. 

Comparing turbulent transport coefficient turbulent resistiv- 
ity with the B z (y) case, we note that no turbulence weakening 
occurs in the present case, consistently with the linear justifica- 
tion (the expected linear growth rate is not modified by a varying 
toroidal field). Although other experiments for larger values of B 
would allow a better determination of the resistivity, we can still 
approximate 



r] zz ~ 0.57a- 

with a mean vertical field and 



'/ 



1.7a 



(38) 



(39) 



0.235a- 



(34) with a mean toroidal field. 
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Fig. 4. Mean field (left) and Emfs (right) from run ZZY4 (B z (y) case) with a mean vertical field B° z = 0.1 and 6B° = 0.04. The right 
plot shows the best fit corresponding to the turbulent resistivity model. We measure in this case r\ xx — 5.6 x 10 -2 . Note that the mean 
field fi? has been subtracted out. 



0.06 



-0.02 





Fig. 5. Mean field (left) and emfs (right) from run ZXY5 (B x (y) case) with a mean vertical field B° = 0.1 and 6B = 0.06. The right 
plot shows the best fit corresponding to the turbulent resistivity model. We find in this case rj a — 8.2 x 10~ 2 . 



4.3. Toroidal field with a vertical dependency B x (z) 
We define the mean field by: 

(B z ) = B° z , 

(B x ) = B° + ( SB°co S (|^). 



(40) 
(41) 



model 


B u 


6B U 




a 


ZXZ1 


0.1 


0.01 


4.8 x 10~ 3 


4.4 x 10~ 2 


ZXZ2 


0.1 


0.02 


2.5 x 10 3 


4.3 x 10~ 2 


ZXZ3 


0.1 


0.03 


3.7 x 10~ 3 


4.6 x 10~ 2 


ZXZ4 


0.1 


0.04 


4.8 x 10~ 3 


4.8 x 10~ 2 


ZXZ5 


0.1 


0.06 


6.8 x 10 3 


5.1 x 10- 2 


ZXZ6 


0.1 


0.08 


3.2 x 10~ 3 


7.2 x 10~ 2 



Table 3. Main results from the B x (z) (toroidal field with a verti- 
cal dependency) case with a mean vertical field 



The mean current is then found to be: 



J = J y e z = d,B x (42) 

(J) = (J y )e y = -SB^j- sin (^) (43) 

Looking for a correlation between J y and £ is more prob- 
lematic in this case. Indeed, if B? + 0, the imposed B x (z) ex- 
cites a channel flow solution, leading to the production of a large 
scale By(z) (see Fig. [6). The existence of a persistent channel 
flow solution also leads to a double period S z (z), as observed in 
the simulation. 

These complications induce the existence of off-diagonal 
components of the resistivity tensor, which cannot be quantified 
with the simple method adopted in this paper. Nevertheless, in 



spite of these difficulties, the procedure used in the previous sub- 
sections yields reasonable results for the diagonal component, as 
can be seen on Fig. [6] 

T]yy ~ 0.08or. (44) 

This value of rj is more uncertain than the other ones derived in 
this work, on the order of 30% to 50%. The correlation with a is 
also not quite satisfied (see next subsection), but remains accept- 
able at the level of precision of determination of this quantity. 

The value of 77,,, is more precisely derived when = 
which prevent channel flow formation. We present on Fig. [7] an 
example of the profiles obtained from such a simulation and the 
resistivity values for several runs in Tab. @] 
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model 




B u 


SB U 


Vyy 


a 


XXZ1 


0.1 


0.0 


0.01 


5.1 x 10~ 3 


6.6 x 10~ 3 


XXZ2 


0.1 


0.0 


0.02 


5.0 x 10 3 


7.8 x 10~ 3 


XXZ3 


0.1 


0.0 


0.03 


5.6 x 10 3 


7.6 x 10~ 3 


XXZ4 


0.1 


0.0 


0.04 


5.5 x 10 3 


8.1 x 10~ 3 


XXZ5 


0.1 


0.0 


0.06 


5.1 x 10~ 3 


9.6 x 10~ 3 


XXZ6 


0.1 


0.0 


0.08 


6.2 x 10~ 3 


1.2 x 10~ 2 


XXZ7 


0.1 


0.0 


0.1 


6.0 x 10~ 3 


1.3 x 10~ 2 



Table 4. Main results from the B x (z) (toroidal field with a verti- 
cal dependency) case for various toroidal B° v 



Interestingly, r] yy seems to be constant although a varies sig- 
nificantly, contrarily to previous cases. We can still approximate: 

rjyy ~ 0.12a (45) 

We note however that this relation doesn't mean that the dif- 
fusive process is very strong in the vertical direction compared 
to that in the horizontal direction, as the turbulence intensity is 
weaker with a mean toroidal field than with a mean vertical field 
(compare for instance the values of a in Tab. [4] with models ZXY 
in Tab.|2]i. However, comparing only mean toroidal field simula- 
tions (X** models), we can state that on average the rj yy coeffi- 
cient is significantly smaller than rj zz . 



4.4. Parameter and configuration dependence of the 
turbulent resistivity tensor: 

It is of interest to compare the effect of the field configuration 
and magnitude on the various turbulent resistivity tensor compo- 
nents that we have characterized, most notably concerning the 
radial and vertical diffusion of vertical and azimuthal field com- 
ponents. To this effect, the various dependencies have been rep- 
resented on Fig. [8] 

The most significant trend is due to the direction of the mean 
field. As for turbulent viscous transport (but somewhat less dra- 
matically), the turbulent resistive transport in the radial direction 
is substantially reduced when the mean field is azimuthal with 
respect to a vertical mean field, by nearly an order of magnitude; 
however, and somewhat surprisingly, no such trend is visible for 
the diffusion in the vertical direction (bottom left subplot). Our 
simulations with a mean toroidal field may suffer (as all others) 
from a lack of resolution of the relevant small scale modes, a 
point that has apparently never been adequately checked in the 
literature. 

Another important trend, and the most critical one for the 
issue of jet driving from accretion disks, concerns the relative 
efficiency in the radial diffusion of vertical and azimuthal field 
components (top right subplot). The diffusion of the azimuthal 
field is more efficient by a factor of order 2 to 3. 

There is apparently some influence of the magnitude on the 
mean field on the efficiency of transport in the radial direction 
(top left subplot). This trend is the same as for the turbulent vis- 
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Fig. 8. Dependence of the turbulent diffusivity tensor components that we were able to measure on B (relative field strength) and e 
(relative amplitude of the varying field component), for various field configuration. The top left figure represents rj xx for two different 
values of B. The bottom left one shows the role of the mean field orientation on tj zz . The top right compares r\ xx and r\ zz . The bottom 
right figure shows the ratio of the relevant turbulent r\ component over the turbulent viscosity. The symbols are the same as on the 
other subplots. 



cous transport, as shown on the bottom right subplot. The sim- 
plest explanation for this trend is that the transport is magneti- 
cally driven, whereas the usual adimensionalization of quantities 
(which we have more or less followed in our definitions) scales 
transport with the sound speed and the scale height: within fac- 
tors of order unity, v = ac s H ~ a£lH 2 and similarly for the 
various rj components. With this scaling, the various simulations 
available in the literature imply that a 3/T 1 ^ 2 . If a scaling with 
the Alfveen speed were chosen instead, i.e., v = (YaVaH, the 
characteristic dimensionless number a a is independent of B (at 
least until Bq is so weak that the zero mean field MRI dynamo 
process takes over). It has been noted in the first subsection that 
the decrease with e for the largest value of B is most likely a sat- 
uration process linked to the fact that this choice of B is close to 
the marginal stability limit of the instability. 



The last subplot (bottom right) indicates that B (amplitude 
of the mean field) and e (amplitude of the field variation) influ- 
ence in much the same way the turbulent resistive and turbulent 
transport. To some extent, this is also true of the spatial orienta- 
tion of the underlying quantities, and not only their magnitudes. 
Overall, the turbulent resistive transport has an efficiency which 
is smaller but comparable to the turbulent viscous one, within a 
factor of order 2 or 3. 



5. discussion 

We have presented a systematic method to determine the tur- 
bulent resistivity associated with MRI turbulence in accretion 
disks. We have exemplified this method in the configuration ra- 
dially and vertically varying vertical magnetic fields, using non- 
linear spectral simulations of turbulence. We have also analyzed 
the resistive transport due to channel modes in the linear limit, 
in the presence of a radially varying vertical fields. Both the lin- 
ear calculations and the fully nonlinear simulations show that the 
turbulent resistive transport is large, comparable to the turbulent 
viscous transport, albeit a factor of 2 to 4 smaller. This feature 
contradicts the heuristic model developed by Shu and coworkers 
dShu et alj|2007l) for turbulent resistive transport in YSOs accre- 
tion disks, thereby removing an obstacle to the launching of jets 
by accretion disks. 

More precisely, we find that the turbulent resistivity rf is 
largely correlated to the turbulent viscosity v T over a wide range 
of variations of the dimensionless parameters of the problem (but 
we have not explored the dependency on the dissipation num- 
bers, which needs an extensive study in itself). As a matter of 
fact, we may define a turbulent Prandtl number Pm T = v T /rj T , 
which is found to be on the order of 2 - 5. Second, we find that 
the turbulent resistivity is an anisotropic tensor. In particular, the 
toroidal field (B x ) diffuses about 3 times more rapidly than the 
poloidal field (Z? ), in the radial as well as the vertical directions. 
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We also find that non diagonal te rms of the turbulent resi stivity 
tensor are non zero. As shown by Lesur & Ogilvie (2008), such 
terms might play an important role for disk dynamos and large 
scale magnetic field generation. 

It is often argued that a turbulent Prandtl number vrlrjT 
on the order of R/H is required for turbule nt disks to be 
able to launch jets (se e van Ballegooijen in lBelvederdfT989l 
and iLubow et al.l [l994h . In fact, because the accretion veloc- 
ity in jet-driving disks is larger than in standard accretion 
disks, this requirement overestimates the necessary turbulent 
resistivity, which turns out to be comparable to the turbulent 
viscosity in self-consis tent acc retion-ejection models (see e.g. 
Casse & FerreiraliOQ O: see also Roths tein & Lovelacel2008l and 



Lovelace et al 



2009 for a simplified version of the same argu- 



ment). This large turbulent Prandtl number argument is also of- 
ten invoked to justify that an outer standard disk cannot trans- 
form into an inner jet-launching one in the accretion process, but 
the validity of this conclusion relies heavily on the assumed ver- 
tical structure of the models considered (in particular the mag- 
netic structure of the corona), a point that has not been appropri- 
ately taken into account in the literature up to now. 

On the basis of these results, it seems quite plausible that 
accretion disks have the ability to launch non stationary jets. 
Although the turbulent resistivity we find is somewhat too weak 
to allow for the existence of stationary accretion-ejection struc- 
ture, the anisotropy is in the right range. Nevertheless, fur- 
ther work is required to get a complete characterization of the 
turbulent resistivity. In particular, the correlation with turbu- 
lent viscous transport needs to be more precisely studied, as 
well as the impact of the (molecular) Prandtl number, which 
is known to be strong on the momentum transport efficiency 
dLesur & Longarettill2007l) . 

While we were writing this paper, a similar study in shearing 
box with the ZEUS cod e has appeared on the astro-ph ArXiV 
by Guan and Gammie dGuan & Gamm ie 2009), and some dis- 
cussion of the connection between the two investigations is in 
order. In their paper, they impose of mean toroidal field, wait for 
turbulence to reach a stationary state. They superimpose then a 
sinusoidally varying field component whose decay rate is used 
to quantify the turbulent resistivity. Instead, we impose a con- 
stant sinusoidal component of the field and study the correlation 
between the resulting emf and the current in the insuing statis- 
tically stationary turbulent regime. Guan and Gammie have also 
made some resolution studies, which we have not performed, as 
previous experience with the dimensionless numbers used in this 
work has shown us that the dissipation scales are adequately re- 
solved with our adopted resolution. Finally, Guan and Gammie 
have looked into the effect of the box aspect ratio, which was 
held fixed here. Conversely, our parameter study is somewhat 
more extensive than theirs for the type of configurations we have 
looked into. The majority of the runs performed by Guan and 
Gammie have been made for vertical sinusoidal component su- 
perimposed on the mean toroidal field, a configuration we have 
not investigated, so that these runs are not directly comparable 
to ours. Note however that using a mean vertical field instead of 
a mean toroidal one is more relevant to the question of flux dif- 
fusion for accretion-ejection m odels. Nevertheless, some of the 
runs of lGuan & Gammiel d2009l) have been performed with both 
the mean and sinusoidal component in the azimuthal direction, 
which are comparable to the ones presented here in section 4.3. 
Their turbulent resistivity and viscosity values are systematically 
larger than ours, but note that their e parameter in these runs is 
four times larger than the larger one we have used (and believe 
to be more relevant to disk physics). As there is a clear trend 



towards a significant increase of the turbulent resistivity and vis- 
cosity with increasing e at the larger values of e shown in table 
4, we conclude that the two results are reasonably compatible. 
More generally speaking, in order of magnitude, both studies 
are rather complementary and agree on the fact that the turbu- 
lent resistivity and viscosity are comparable, but runs in similar 
settings would be required to make a full comparison between 
the two methods. 
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